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THE SLIDES (C18) 

The four slides in this set were all created by Ken Philip. 

457: Center = —1.71115. 

535: SHVRBS58.fs. 
Center = —.7459,0601,58 + .0984,8530,361i, magnifica- 
tion = 10°, escape radius = 2, dwell limit = 3825, program 
= MandelZot 1.3. A complex ‘quad spiral’ structure 
around a midget deep it Seahorse Valley. The previous 
batch had other views of this complex object. The effects 
of aliasing can produce surprising changes in the per- 
ceived patterns as the color mapping is altered, so it is 
worth experimenting with color mapping here. 

543: First Pinwheel. 
Center = —.1986,1627,510 + 1.1003,4606,584i, magnifica- 
tion = 5x108, escape radius = 2, dwell limit = 1000, pro- 
gram = MandelZot 1.3. A rerun in MandelZot of the first 
pinwheel I ever found, in a galaxy cluster on a tendril on 
the midget above the N wart. #543 uses a smooth CLUT, 
with different settings of the color animation cycle. 

544: Furry galaxies. | 
Center = —1.7219,0586,9244 659 — .0022,0479,8027,7931, 
magnification = 10", escape radius = 2, dwell limit = 
1000, program = MandelZot 1.3. A galaxy cluster off a 
midget on the spike, decorated with asymmetrically dou- 
bling spines. 


SLIDE #547 
Slide #547 was distributed with slide set C17, but only the 

name was given in Amy #17. Here is some more informa- 

tion about it. 

547 (Ken Philip): Furry galaxy pinwheel.l. Center = 
~1.7219,0586,9244 66016 — .0022,0479,8027,794761, 
magnification = 8.7x10’* escape radius = 2, dwell limit = 
1000, program = MandelZot 1.3. A detail of the center of 
a galaxy cluster (#544) off a midget on the spike, decorat- 
ed with asymmetrically doubling spines. 
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MORE SZILAGYI SLIDES 


Steven Szilagyi- -has kindly supplied two more zoom se- 
quences of 20 color slides each. These are available from 
Amygdala: see the Price List / Order Form for prices. 


RENEWAL TIME! 

Look at the label on the envelope this newsletter came in. 
In the lower right there is a box with one or two numbers in 
it, indicating when your newsletter subscription expires and 
when your slide subscription expires (if you get it). N20, for 
example, means that newsletter #20 is the last issue you will 
receive per your current subscription. 

C16 means that slide set #16 was the last one you re- 
ceived. If that is what you’ve got, it means that as of this is- 
sue (#18) you’ve already missed two slide sets, #17 and #18! 
If you want to renew, here are the prices (postpaid) AS OF 
THIS DATE. Slide sets 17-22 (or any 6 sets, 4 slides each): 
$19.20; or slide sets 17-26 (or any 10 sets, 4 slides each): 
$32.00. If you’ve reading this much after September, 1989, 
please refer to a current price list. 
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FRACTALS AND HDTV 

A short piece entitled GEOMETRY COULD GIVE HDTV 
SIGNALS THE RIGHT SHAPE appeared in Business Week 
for April 3, 1989, page 110: 

The wide-screen, supercrisp pictures promised by high- 
definition TV require much more information than ordinary 
telecasts carry. Trying to squeeze an HDTV signal into exist- 
ing broadcast channels, as the FCC has mandated, is giving 
many researchers grief. But mathematicians Michael F. 
Barnsley and Alan D. Sloan at Georgia Institute of Technolo- 
gy may have a solution: fractal geometry. 

Fractals are those seemingly haphazard lines that turn out 
to have a strong kind of order. They are helping scientists 
gain insights into patterns hidden within meandering shapes 
such as coastlines and mountaintops. Barnsley and Sloan say 
the esoteric formulas can be used to encode an image — and 
squeeze it to a mere 1/1000 of its original size. 

They believe the approach offers a way to compress a TV 
signal for transmission, then reverse the process in the TV 
set. Barnsley hopes to develop a system by year end. Mean- 
while, the duo’s Atlanta-based startup, Iterated Systems Inc., 
in February started marketing a $32,500 software package for 
creating models of biological systems from fractals. 

— Sent in by Richard Morley 


High-definition pictures without HDTV signals?, from Ma- 
chine Design(May 25, 1989) p 15. 

An Atlanta-based company called Iterated Systems Inc. 
hopes to generate high-definition television (HDTV) pictures 
from an ordinary broadcast signal. The trick is an automatic 
enhancement system which generates high-resolution images 
in real time. ISI’s system will use a black box installed on 
HDTV receivers. Since the new technology is only on the re- 
ceiving end, the usual problems associated with HDTV 
broadcasting are avoided. Says ISI President Alan Sloan, 
“We require no changes to the transmitted signal, so our tech- 
nique is compatible with current NTSC broadcasting. More- 
over, Our system can be used with any of the HDTVs current- 
ly being proposed, so we are not worried about the adoption 
of standards.” 

Iterated Systems’ technique is based on fractal geometry, 
which converts the incoming signal to a fractal representation 
in the form of a digital code. The fractal representation is in- 
dependent of the display device and so the image can be 
viewed as video having higher resolution. Says Michael 
Barnsley, executive vice-president of Iterated Systems, “Be- 
cause the fractal representation has no scale attached to it, 
you can recreate the image to whatever resolution you want.” 

A fractal is basically a complex geometric shape with 
very low information content. “In other words,” says Barn- 
sley, “fractals can produce visually complex and realistic im- 
ages and yet can be specified and controlled with minimal 
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amounts of data.” 

One result of Iterated Systems’ work toward high- 
definition television is a Video Rate Decoder that is expected 
to debut this summer. Iterated Systems intends to store 45 
seconds of video onto a high-density floppy disk that can, 
through the decoder, be played back on a PC/AT computer. 
Says Barnsley, “The image is high resolution and the quality 
is good. It’s not HDTV but it does use a similar technology.” 
— Sent in by A. R. Wadum 


DWELL FREQUENCY PEAKS IN MIDGETS 
ON THE SPIKE 
— Ken Philip 

In the description of color slide #144 (by Andrew La- 
Mance) in AMYGDALA #13, there was a diagram showing 
peaks (for every fourth dwell value) in a dwell frequency plot 
for the midget at -1.94254 on the spike (See Figure 1). This 
behavior struck me as interesting, and I decided to look into 
it further. [Note: The center of the midget itself is approxi- 
mately at -1.942525. -1.94254 refers to a point a bit to the 
left of the center of the head of the midget.] 


The first test I made, after modifying a B&W True BASIC 
Mandelbrot program to produce dwell frequency histograms, 
was to verify LaMance’s plot for this midget when the es- 
cape radius was set at 2, and then replot the histogram with 


LaMance’s Midget 


Escape radius = 2 


— Peaks every 4th dwell step 


Aj 


Escape radius = 100 


Figure 1 — Dwell Frequency Plot 
Re=-1.94254, Im=0, mag=13221, max dwell=256 
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an escape radius of 100. I suspected that the peaks might be 
an artifact of using 2 as the escape radius, and thus was not 
surprised to see no trace of the peaks on the ER=100 plot 
(see Figure 1 on the previous page for a comparison of the 
two histograms). I then ran two B&W prints of the midget 
with escape radii of 2 and 100 (Fig. 2) to see what on the im- 
ages corresponded to the dwell peaks. It was clear (Fig. 2a) 
that the bands in the ER=2 image with the very wide lobes 
were responsible for the dwell peaks in the histogram. 


LaMance’s Midget 


- TOATE = 
14 vA 
slater 


Re=-1.94254, Im=0, mag=13221, dwell=256 
Escape radius=100 


LaMance’s Midget 


Figure 2a 
Re=-1.94254, Im=0, mag=13221, dwell=256 
Escape radius=2 


Amygdala #18 


Figure 3 
Spike: showing dwell bands 
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LaMance midget 


John Dewey Jones midget 


Figure 4 
Spike tip: showing dwell bands 


Figures 3 and 4 (showing the entire spike, and a magnified 
view of the end of the spike) indicate that the LaMance mid- 
get is formed by a compression of the 4'™ dwell contour (in 
an ER=2 plot). It occurred to me to wonder if other midgets 
would manifest the same behavior, and if there was any cor- 
relation between the number of the outermost dwell contour 
that was compressed around the midget, and the spacing of 
the dwell frequency peaks. Figure 5 (next page) shows the 
dwell frequency plot (ER=2 and ER=100) for the John 
Dewey Jones midget, formed by a compression in the 6" 
dwell contour, and Figure 6 (next page) shows the actual 
midget with ER=2 (Fig. 6a) and ER=100 (Fig. 6b). It is clear 
that the dwell peaks occur for ER=2, but at every 6" dwell 
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John Dewey Jones’ Midget 


John Dewey Jones’ Midget 


Escape radius = 2 


every 6th dwell sep 


dwell 100 


Escape radius = 100 


Figure 6b 
Re=-1.99638, Im=0, mag=150000, dwell=256 
Escape radius=100 


sion in the 5“ dwell contour. Figure 7 shows the results: the 


Figure 5 
Dwell Frequency Plot peaks are not so pronounced as with the two preceding mid- 
Re=-1.99638, Im=0, mag=150000, max dwell=256 gets, but they are present, and at a spacing which correlates 


with the outermost dwell band compression. 


Biggest Midget - 1.861534 midget 
Escape radius=2 Escape radius=2 


Peaks every 3rd dwell sep Peaks every Sth dwell sep 


Figure 6a 
Re=-1.99638, Im=0, mag=150000, dwell=256 
Escape radius=2 


Figure 7 
Dwell Frequency Plots 


value. Again, for ER=100 there are no such peaks. 
Figure 8 (next page) shows the more pronounced peaks on 


I ran dwell frequency plots (ER=2) on the Biggest Midget, | the midget at -1.985444, with a peak at every 5" dwell con- 
which is formed by a compression in the 3% dwell contour— | tour. This midget is also formed by a compression of the 5" 
and on a midget at -1.861534 which is formed by acompres- | dwell band. 
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Biggest Midget 


- 1.861534 midget 


Escape radius=2 Escape radius=2 


Peaks every 3rd dwell step Peaks every Sth dwell stp 


Figure 7 
Dwell Frequency Plots 


-1.985444 midget 


Escape radius = 2 


Peaks every Sth dwell sep 


Figure 8a 
Dwell Frequency Plot 


The cycle number of the Biggest Midget is 3, and a run 
with a modified version of my WormTracks program indicat- 
ed that the John Dewey Jones midget may have a cycle num- 
ber of 6. Before one jumps to the conclusion that the cycle 
number is always equal to the dwell number of the outermost 
band compression, and in turn controls the dwell frequency 
peaks, I should mention that a WormTracks run showed 5, 
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rather than 4, fixed points for the LaMance midget (and one 
of these was rather wide, suggesting that under higher magni- 
fication it might be double, which would give a total of 6 
fixed points). So the relationship between the cycle number 
and the dwell peaks does not seem quite clear yet, although 
there may be such a correlation for midgets at or near the 
ends of their ‘series’ (as: the JDJ midget is the last midget 
formed by a compression of the 6" band). 


I leave it to the mathematicians to explain why an escape 
radius of 2 gives rise to this phenomenon... It is clear, how- 
ever, that using an escape radius of 2 means that you are re- 
stricting the entire ‘universe’ to the inside of a circle about 
the origin with that radius — since you can obtain no data 
from beyond the escape radius with a conventional Mandel- 
brot program. On a purely intuitive basis, using an escape ra- 
dius large compared to the size of the Mandelbrot Set would 
appear more useful, since that allows you to follow dwell 
contours out to large distances. 


-1.985444 midget 


Figure 8b 
Re=-1.985444, Im=0, mag=6000, dwell=200 
Escape radius=2 


TUTORIAL — THE BINOMIAL EQUATION 
(PART 2) 
— RS 
(Continued from Amy 16) 

In this installment we are going to give solutions to the ex- 
ercises posed at the end of the last installment. 
1. Express cos 40 and cos 5@ in terms of cos 8 and sin 8. 

De Moivre’s formula can be used to express cos 48 in 
terms of cos @and sin 8. Taking n = 4, (16.3) becomes: 


cos 40 + i sin 480 = (cos 0 + i sin 6)* (18.1) 
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The binomial formula (16.4) applied to the case n = 4 
gives: 


(a+ b)*= > (4 )a"b*” 


The coefficients (4). pe o) are 1, 4, 6, 4, 1, so (18.2) be- 
comes: 
(a+b)* = a‘ + 4a°b + 6a°b’ + 4ab° +b’ (18.3) 
Using (18.3) to expand the right hand side of (18.1) we 
get: 


(18.2) 


cos 40 + i sin 40 = 
cos*@ + 4i cos?@ sin 8 -— 6 cos?8 sin2@ 
— 4i cos 8 sin?0 + sint8 (18.4) 
If we equate the real parts of the two sides of (18.4) separ- 
ately, we get: 
cos 48 = cos*@-— 6 cos2@ sin?0 + sin*0. 


The case n = 5 can be treated similarly. The binomial co- 
efficients are 1, 5, 10, 10, 5, 1, and: 
cos 50 = cos°@— 10 cos?@ sin2@ + 5 cos @ sin*@ 


The following identity will prove useful in solving prob- 
lems 2 - 5: 

(2+... 2002) = 1-2" (18.5) 
2. Simplify 1 + cos O + cos 20 + ... + cos nO and 
sin Q + sin 28 +... + sin nO. 

If 8 is a multiple of 27, we have 1 + cos 6+ cos 20 + ... + 
cos nâ = n+1, and sin 6+ sin 20 + ... + sin nð = 0. Let us 
next suppose that @ isn’t a multiple of 27. 

Using (16.3) for n = 2, 3, ..., we get: 

1+ (cos 0 + i sin O) + (cos 8 + i sin 6)? +... + 
(cos 0 + i sin 6)” = 
1 + cos 0 + cos 20 + ... + cos nO + 
i(sin 8 + sin 20 + ... + sin n8) (18.6) 

Let z = cos 0 + i sin 0, and note that z + 1, since 8 is nota 

multiple of 27. Thus from (18.5), 
(l+¢zt+..t¢2") = (1-2) -2) 


or 
1+ (cos 8 + i sin 6) +...+ (cos O + i sin 6)” = 
(1 — (cos 0 + i sin 6)"*')/(1 — cos @ — i sin 8) 
(18.7) 
(18.6) and (18.7) yield: 


1 + cos 0 + cos 20 + ... + cos nO + 
i(sin 0 + sin 20 + ... + sin nO) = 
(1 -= (cos 8 + i sin 6)"*')/(1 — cos 0 — i sin 6) (18.8) 
Equating real and imaginary parts of (18.8): 

1 + cos 0 + cos 20 + ... + cos n = 

R[(1 - (cos 8 + i sin O”**)/(1 — cos 0 — i sin 6)]; 
sin 8 + sin 20 + ... + sin nô = 

3[(1 —(cos 8 + i sin 6)"*!)/(1 — cos 0 — i sin 0)]. 

Does that qualify as a “simplification”? 
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3. Express the fifth and tenth roots of unity in algebraic form. 
Since 1 is itself “a” fifth or tenth root of unity, I assume 
we’re being asked to find primitive roots of unity. 


Fifth root: @s5= cosaz +i sin2Z. 


Using (16.3) and letting c =cos22, s = sin2Z: 
5 „5 5 5 
l= œo = (c +is) = 
© + 5ic4s — 10c3s* — 10ic?s? + Scs* + is? 

so 0 = 5c*s— 10c?s? + s’. We can divide by s, since 
s#0: 0 = Sc*- 10c*s* + st. But s? + c? = 1, so 
0 = 5(1-s”)? — 10(1-s”)s? + s* = 1654 — 20s? + 5. The result 
1S: 


(18.9) 


16s* — 20s* + 5 = 0 (18.10) 
Whew! What started out looking like an intractable fifth 
degree equation, then a monstrous fourth degree equation, 
has tumed into a pussycat quadratic in s*. Solving (18.10) by 


2 SHS ors =t pace 


the usual method, s^ = Which 


root? 
Well, consulting a trig table or a calculator, s = sin2% = 
sin 72° is about 0.95, so we want the positive root: s = 


5+¥5__ We then have c = Vl-s* = +4 | 1 - 835, = 
8 


tA | s , and since c = cos 72° = 0.31, we take the posi- 
tive root: c = vero 


Finally, œ, = c + is, which in algebraic terms is: 


Tenth root: As far as @ = cost + isin2® goes, the same 


approach will not work, since the equation analogous to 
(18.9) is of too high degree to solve. Fortunately we have, by 
(16.3), @,* = coset + isin? = @;, SO We get: 


8 8 


We need to tidy this up, i.e. to express @,, in the form a + 


ib. If œ =a + ib, then (œ) = (a + ib)’ = a° — b? + 2iab. 


But @ = -4/28 vin [S215 so we have: 


(18.11) 


a? — b? = 
ab Seppo. =<. (18.12) 


(18.12) gives us b = oy. which substituted in (18.11) 
a 
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A 2 
gives: A = c. Thereforea* — ca? ve = 0, whose solu- 


2a 
tion is a? = &+L, so a =4/e+l =4/ 14 3-V5 and 
2 2 2 32 
b=S =4/1- oD... Summing this up: 
2a Z 32 
Hp: = Lo 3-¥5 +i L 3-15 
32 Z 32 


4. Ifo= cos2Z + i sin2Z, prove that 


1+ ot + oht.. ol™™ = O, if his an integer not a mul- 
tiple of n. 
Let z = œ. Note that: (1) œ” = 1, and (2) since A is an inte- 


27h 


ger not a multiple of n, cos=" cannot be 1. Applying 


(16.3), z = (cos22 + i sin22)' = cosź% + i sin 2 oe 
Therefore 
1+ o+ wo +...+ ol) = 
1+z+z?+... +z") = 


n hn nh h 
i E Ore ER a TE, ak ig 


5. Whai is the value of 


Ss 1-@"+ o*-...+ Aye? 
Letting z = —a, the sum becomes 
= l+z¢7+...¢2"7 (18.13) 
If af =—1, i.e. if z= 1, then S is the sum of n ones, i.e. 
if af = -1 then S =n. (18.14) 


Let’s deal with the other case: œ” + —1, i.e. z #1. Recall- 
ing that w is an n-th root of unity, i.e. that œ” = 1, we have 
z” = (-w")” = (-1)" (@")" = (-1)".. Applying (18.5) to (18.13): 


z"-1 _ 1-(-1} 


S =“ = ~~ soif nis even, S = 0, while if n is 
odd, 2 FA u Combining these results with (18.14) we 
can sum up: 


If œ = —1 then S = n; 
if wo + —1, then if n is even, S = 0, 


Z 
while if n is odd, 9 = -3 
w+ 1 
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DIVCON ENHANCEMENTS 
— Dave Platt 

I've made a number of enhancements [in MandelZot] to 
the Divide and Conquer code. I first added some calls to the 
point-drawing code, so that the user can observe DivCon’s 
progress around the perimeter of the rectangle under exami- 
nation. Once I saw what DivCon was actually doing, I real- 
ized that a simple and effective optimization was possible, 
and incorporated the optimization into the code. 

Briefly: before the code scans along one of the sides of a 
rectangle and computes the dwells, it performs a quick 
glance along the next “outermost” row or column. If all of 
the points immediately “outside” of the rectangle’s side have 
already been calculated, and if their dwells match the dwell 
of the rectangle’s corners, then the code skips the calculate- 
and-check pass on that side, and goes on to the next side. For 
example: 


----AAAAAAAAAAAAAAAAAAAAAAAA---- 


--LS------------------------ SR-- 
--LS------------------------ SR-- 
--LS------------------------ SR-- 
--LS------------------------ SR-- 


---CSSSSSSSSSSSSSSSSSSSSSSSSC--- 
----BBBBBBBBBBBBBBBBBBBBBBBB---~- 

In this case, assume that the rectangle being scanned is de- 
noted by the S (side) and C (corner) characters. The algo- 
rithm checks the four C corners to see if their dwells are the 
same and exits (or recurses) if not. If the C’s are the same, 
then the algorithm quickly scans the A (above) pixels; if their 
dwells have all been calculated and are the same as the C 
dwells, then the program assumes [provisionally] that all of 
the top-side S pixels would also have the same dwell, and 
doesn’t bother to calculate them. If any A pixel hasn’t been 
calculated, or doesn’t match the C pixels, then the program 
runs along the top-side S pixels, calculating each one (if not 
already calculated) and comparing with C. An analogous 
check is made on the right (R), bottom (B), and left (L) sides. 

This modification is a real win when the program is Div- 
Con’ing an area that includes a large chunk of the M-set. If 
the rectangle being scanned lies entirely within M, then the 
program won’t need to calculate the top-edge or left-edge 
pixels for this rectangle — the ones adjacent will already 
have been calculated in an earlier phase of DivCon. Net time 
savings for the rectangle in question: 25% to 50%. 

In a “vanilla” DivCon run (in which the entire window is 
being scanned, and no points in the window have been 
scanned before), at most two sides of any rectangle-under- 
test will be adjacent to other already-scanned rectangles. In 
such runs, there’s no advantage to performing the adjacent- 
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row/column scan for the other two sides (bottom and right, in 
MandelZot’s case). However, some programs may permit 
the user to select a portion of the window and run DivCon on 
only that portion. If the user does this several times to adja- 
cent or overlapping portions of the window, then it’s possible 
for a rectangle-under-examination to have three or (in degen- 
erate cases) four sides adjacent to already-calculated portions 
of the window. 

The scan of the row/column adjacent to the rectangle un- 
der consideration is very quick, because the program makes 
no attempt to calculate dwells for the points being scanned. 
This reduces the cost of a “failed” scan to a very small frac- 
tion of the benefit of a “succeeded” scan. 

I’ve made one other small optimization to DivCon. Be- 
cause the M-set and its surroundings are symmetrical with re- 
spect to the Y=0 axis, there’s no need to run DivCon on pix- 
el-rows where Y<0 if it has already been run for the 
corresponding Y20 rows. The top-level DivCon code in 
MandelZot now runs DivCon for the rows that cover Y20, 
and then copies the data to the corresponding rows for which 
Y<0. This speeds up the DivCon quite a bit when examining 
the whole M-set, or when looking at structures that straddle 
the axis of symmetry (e.g. the John Dewey Jones midget). 

Dave Platt 

1095 Burgoyne St. 

Mountain View, CA 94043 
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Springer. I [Reinald Eis] read a couple of these volumes and 
would like to recommend a few articles to the readers of 
Amygdala, because the topics dealt with pertinent fractals, 
too. I ought to include a cautionary word to readers who are 
easily frightened off by maths, even though most articles fea- 
ture nice illustrations to improve digestibility.] 


PRODUCTS 


1990 FRACTAL COSMOS CALENDAR 

The 1990 Fractal Cosmos calendar, published by Am- 
ber Lotus, is now available from Amygdala. It features 
13 beautiful fractal images in color done by Peitgen and 
colleagues. 

The price is $9.95 plus $2.50 postage in the U.S. (First 
Class), Canada, and Mexico (Air Mail). Air Mail post- 
age elsewhere is $4.90. 


J.E. LOYLESS: ; 5185 Ashford Ct. / Lilburn, GA 30247 
1. DISCOVER CROSS-EYED 3-D 

If you can cross your eyes, you can see and manipu- 
late true 3-dimensional split-image wire models of poly- 
hedra, and nest them in the manner of M.C. Escher. For 
IBM compatibles w/ EGA or Hercules. BASIC included. 
$5.00 postpaid. 


2. NEW FRACTALWARE 
These fractal routines are written for flexibility rather 
than speed, and allow a wide variation on the Mandelbrot 
theme, including true 3-D, perspective, and interior 
views using four formulae in any of six possible planes. 
For IBM compatibles w/ EGA or Hercules. BASIC in- 
cluded. $6.00 postpaid. 


3. 1990 FRACTAL CALENDAR 

Thirteen novel laser-printed monochrome fractal im- 
ages, including perspective and interior maps of the Man- 
delbrot set. $5.00 postpaid. 


CIRCULATION 

As of September 24, 1989, Amygdala has 809 paid-up 
subscribers, 178 of whom have the supplemental color slide 
subscription. 
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